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ABSTRACT 


High speed oscillatory flows, like high speed mean flows, are capable of inducing 
time-averaged heat transfer effects. This research involves the analytical solution of a model 
problem of zero-mean internal oscillatory flow, which arises from a high-intensity resonant 
standing acoustic wave set up across the ends of two parallel plates. The compressible form 
of the Navier-Stokes equations are solved, along with the equations of continuity, energy, 
and state, using perturbation solution and complex variable methods. MAPLE, a symbolic 
mathematical software tool, is utilized to find the time-averaged portion of the temperature 
distribution between the plates. The final heat transfer results are presented in terms of 
suitably defined recovery factors. The analysis is performed for varying gap widths between 
the plates using air as the host fluid. This work provides the fundamental explanation of the 
phenomenon responsible for the thermoacoustic refrigerating effect as well as an analytical 
basis for determining the optimum gap width between the plates of the stack in a 
thermoacoustic refrigerator. 
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I. INTRODUCTION 


This research involves the study of the time-averaged heat transfer effects induced 
by high speed internal zero-mean oscillatory flows. The model problem being treated is 
for an ideal-gas host fluid between two “infinite”, parallel plates. The zero-mean 
oscillatory flow is produced by a high-intensity standing axial acoustic wave induced 
across the ends of the plates. The plates are separated by a distance, 2y 0 , and the sound 
wave has a wavenumber of k=2n/X where X is the wavelength. The relationship between 
these two parameters, y 0 and k, is assumed to be such that yo&«l throughout the analysis. 
In other words, the narrow acoustic waveguide approximation, which is discussed further 
in Chapter m, holds true for this research. Another general assumption made for this 
analysis is that the host fluid between the plates is a gas which is treated as ideal. 

The particular focus of this study is to solve analytically the Navier-Stokes 
equations for this problem. The governing equations are set up and simplified using the 
particular geometry and assumptions pertinent to this particular case. Then, the variables 
are scaled in a standard fashion to reduce the equations of continuity, momentum, energy 
and state into dimensionless forms. An order of magnitude analysis is performed to 
further simplify the equations. Finally, perturbation methods are employed, and like- 
order terms are grouped to form new equations. The equations of most significant order 
are solved analytically using complex variable methods for their respective unknowns. 
The final goal is to determine the resulting time-averaged temperature distribution and 
study its affects on the heat transfer behavior at the walls. MAPLE, a symbolic 
computational software package (see Blachman and Mossinghoff (1994)), is utilized to 
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perform the extensive mathematical manipulations that are required in this analysis. The 
results of the time-averaged thermal behavior that is ultimately of interest are presented, 
by convention, in terms of suitably defined recovery factors. 

The entire analysis is performed for two separate cases, which differ based on the 
width of separation between the plates, 2y 0 , in comparison with a length scale called the 
Stokes’s boundary layer thickness, 8y The Stokes’s boundary layer thickness, also called 
the viscous penetration depth, is defined in Schlichting (1979) in the discussion of 
Stokes’s second problem concerning the flow near an oscillating plate as the thickness of 
the layer of fluid which is carried by the plate. If the quantity, y^, is greater than or 
equal to a certain amount (which is determined in this research and discussed in Chapter 
IV), the spacing between the plates is considered to be a “wide gap.” The first and 
simpler case is analyzed using a wide gap geometry and its respective assumptions. In 
the second case, an arbitrary gap width is considered which lifts some simplifying 
assumptions so that it includes “narrow gaps.” The arbitrary gap case is useful for 
analysis of the time-averaged transport behavior in the stack region of a thermoacoustic 
refrigerator. 
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II. BACKGROUND 


Since the time-averaged heat transfer results of this work are presented in terms of 
conventionally defined recovery factors, it is useful to review the concept of a recovery 
factor. In order to do so, the application of recovery factors to two well studied flows, 
namely Couette and Poiseuille flows, will be presented here as adapted from Burmeister 
(1983). 

Couette flow, as illustrated in Figure 2.1, occurs when fluid exists between two 
infinite parallel plates. 



Figure 2.1 Couette Flow. After Burmeister (1983). 


One plate (the upper plate here) is moving with a velocity Vo, and the other 
(lower) plate is stationary. The plates are at different temperatures, T 0 and T x . If both 
plates were stationary, the temperature distribution between them would simply be linear; 
however, the fluid motion caused by the motion of the upper plate affects the temperature 
distribution. This effect comes from viscous dissipation and is especially significant in 
high speed flows. Applying the following assumptions and boundary conditions to the 
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governing equations of continuity, momentum, and energy, the velocity and temperature 
distributions are determined. 

The assumptions being made are: (1) steady conditions, (2) laminar flow, (3) 

constant properties, (4) no pressure gradient in the x direction, (5) no edge effects or end 

effects, ( 6 ) Newtonian fluid, and (7) no buoyancy or other body forces. 

Boundary Conditions: At y=0: u=0=w (no slip at wall) 

v =0 (no penetration at wall) 

T=T 0 (fixed wall temperature) 

At y=2y 0 : m=V 0 , w= 0 (no slip at wall) 

v =0 (no penetration at wall) 

T=T X (fixed wall temperature) 

The resulting velocity distribution is linear as given below: 


u_ y 
V 0 2y 0 

The resulting temperature distribution is given by 


r-r 0 

T t -T 0 



+ (£cPr)(—^-)( 1 - 

2 y 0 



( 2 . 1 ) 


( 2 . 2 ) 


V 2 

where Ec = ———7 -, which represents a comparison of viscous dissipation effects to 

2 c p (T x -T 0 ) 


the imposed temperature difference, and Pr = —— is the Prandtl number of the fluid. 

k. 


Using Fourier’s law, the heat flux at the lower plate is given by 


PrV 2 

k c [T 0 -{T,+ ——)] 

. dT , i;-r ( rn ^ 2c 

<l w =-k c — = -k c - 7 —-(1 + FcPr) =----— 

d y y=o 2 >'o 2y 0 
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where T waI ^T 0 from before. Now it is observed that the heat transfer is driven solely by 
the temperature difference, T^-T^,, where T m . is not constant but varies depending on the 
particular geometry and flow conditions being studied. 


This result incorporates, into a simple form, the effects of viscous dissipation on 
heat transfer which become more and more significant as V 0 and the kinetic energy of the 
fluid increase. It also allows the heat transfer coefficient to be expressed in a more 
consistent and convenient manner as 
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(T^-T^) 2 y 0 

The heat transfer coefficient, as a result, is never negative; and the Nusselt number, 

2/i 

Nu =--, conveniently has the value of one. 

K 

From the expression for the adiabatic wall temperature in Equation 2.5, it is 

V 2 

clearly observed that some of the fluid’s kinetic energy, —, is being converted into 

2 

thermal energy. In a sense, some of this kinetic energy is being “recovered” as thermal 
energy at the lower adiabatic wall which then results in a rise in temperature of that wall. 
The ratio of the thermal energy recovered at an adiabatic surface to the kinetic energy of 
the fluid is called the recovery factor, r. For this Couette flow, the recovery factor is 

c ,V m -T x ) _ 


The physical significance of the recovery factor, as a measure of viscous 
dissipation effects, applies equally well for other types of flow and wall conditions, even 
for situations with non-adiabatic surfaces. For any situation, “r” contains information 
about how much of the fluid’s kinetic energy is converted into thermal energy at the 
walls. It is defined and calculated, however, in a similar manner to that of the Couette 
flow example - the ratio of the recovered thermal energy at an adiabatic surface to the 
fluid’s kinetic energy. The expression of the kinetic energy may differ depending on the 


velocity chosen to base it on, but the temperature difference used for calculation is always 
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Taw-T wa i h Of course, for other types of flow, the resulting expression for the recovery 
factor may be more complicated than simply a linear Prandtl number dependence. 

Another type of flow is Poiseuille flow which is of direct relevance to the current 
work. In Poiseuille flow, pictured in Figure 2.2, both surfaces are stationary. The 
assumptions applied to this case are the same as for the Couette flow case with the 
exception that there is a uniform pressure gradient in the x direction which drives a uni¬ 
directional mean flow. Using the following boundary conditions, the velocity and 
temperature distributions are determined. 



Figure 2.2 Poiseuille Flow. After Burmeister (1983). 


Boundary Conditions: 


At y=-y 0 : 

u=0=w 

(no slip at wall) 


v =0 

(no penetration at wall) 


T=T 0 

(fixed wall temperature) 

At y=+y 0 : 

u=0=w 

(no slip at wall) 


v =0 

(no penetration at wall) 


r=r, 

(fixed wall temperature) 


The resulting velocity distribution is given by 


i-4) 


2 |X dx 


y o 


(2.9) 
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The definition of average velocity, 


u «v -J-r ud y ’ 

2 y 0 


( 2 . 10 ) 


gives 


3 { dx j 


( 2 . 11 ) 


Also the maximum velocity, which occurs at the centerline (y= 0 ) is given by 


, _ y<> , dp 3 

m 2 ( dx ] 2 Uav 


( 2 . 12 ) 


The velocity distribution from Equation 2.9 may then be expressed in terms of either the 
average or the maximum velocity as 


r 3 v 2 

M = M mO- Z r) or u = - h w (1~) 

2 yl 


(2.13) 


The resulting temperature distribution for this Poiseuille flow is given by 


T=T 0 +^-^(l + ^- ) + ^(l-4) 

2 y 0 3k c y 0 


(2.14) 


From Fourier’s law, the wall heat flux at the lower wall is given by 


k,lT„-(T, +^ 7 %] 

= -K ~r(y = ->o) =--——— 

dy 2y 0 


(2.15) 


Using the concepts established in the Couette flow case, the adiabatic wall temperature 
for Poiseuille flow is deduced to be 


T uw ~T\ + 


6 Prut 


(2.16) 
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This can be also easily verified by solving the adiabatic wall case and finding the 
temperature of the adiabatic wall. 

Again, it is observed that the driving potential for heat transfer between the walls 
is the difference in temperature, T wal rT m . This ensures that viscous dissipation effects are 

k 

accounted for and that the heat transfer coefficient, h = — s ~, is always positive. The 

2^0 


Nusselt number, Nu = , is again unity. The recovery factor, r, is then given by 


, _ p 


c n (T a ~T t ) 


= 12 Pr 


u 


av 

2 


(2.17) 


Note that the kinetic energy of the fluid was chosen to be represented in terms of the 
average velocity. 

These concepts of an adiabatic wall temperature and a corresponding recovery 
factor apply to other geometries and flow conditions and are typically of consequence in 
high speed flows where velocity gradients and viscous dissipation effects are significant. 
In high speed mean flows, a sufficient amount of the fluid’s kinetic energy is transformed 
into thermal energy, causing significant heat transfer to occur at the walls from viscous 
dissipation. Normally, the result is a heating of the walls. 

Recovery factors will be found for the geometry and flow conditions present in 
this research which are similar to the example case of Poiseuille flow just considered. 
The physical geometry is exactly the same with the presence of two parallel plates. The 
principal difference in the current problem is in the flow conditions. In this case a driving 
(zero-mean) oscillatory pressure gradient is set up across the plates by a high-intensity 
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resonant standing acoustic wave. This contrasts to the uni-directional mean flow which 
occurs in Poiseuille flow. Despite being a zero-mean oscillatory flow, the cross-coupling 
between the harmonics is capable of inducing a time-averaged flow and a thermal effect 
which can be significant in high speed situations. This key idea introduces a number of 
interesting new features into the problem which are explored in this research. The posed 
problem and its solution will be addressed further in Chapter HI, and the corresponding 
recovery factor will be discussed extensively in Chapter IV. 
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III. THEORY 


A. THE MODEL PROBLEM AND GOVERNING EQUATIONS 

The model problem consists of high speed oscillatory flow between two “infinite” 
parallel plates and is depicted in Figure 3.1. Both plates are held at a fixed temperature, 
T m . The host fluid, air, oscillates between the plates in the x* direction with a velocity 
amplitude of oscillation, V 0 , and a frequency, to, in radians. 



Figure 3.1 Model Problem: Oscillatory Flow Between Parallel Plates. 


This flow is induced by a high intensity standing acoustic wave set up across the 
ends of the plates. The acoustic wave number is k=2it/X where X is the acoustic 
wavelength. The wavelengths of concern in this research are on the order of 1 meter. 
The distance between the plates is 2 y 0 , and a Cartesian coordinate system is established 
with the origin at the centerline. Using a wavelength of the above mentioned order, 1 
meter, the narrow acoustic waveguide approximation, yoik«l, requires that y 0 be much 
less than 1 meter. In this research, the gap widths (2y 0 ) studied are on the order of 
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millimeters or tenths of millimeters. Therefore, the narrow acoustic waveguide 
approximation certainly holds throughout this work. 


To determine the temperature distribution between the plates, which is the 
fundamental objective of this study, the compressible form of the equations of continuity, 
motion, and energy are solved analytically. These equations in Cartesian coordinates are 
(asterisks (*) denote a dimensional quantity): 

Continuity: —-(p * u *) + -^-(p * v*) + ^(p * w*) = 0 (3.1) 

at * dx * dy* dz* 

Motion: 


x * component: 

9 «* 9 m * 9 m * 9 m * dp* 

p *(—— + u*~ —+ v*-— + w* —— ) = p * gx *—— 
dt * dx* dy* dz dx * 

3 9m* 2 9m* dv* 3w* 

+ a^ B,(2 a^-3 ( a? + ^ + aF )] 

d 9m* dv* d r ,dw* du * x . 


(3.2a) 


y * component: 

n *,dv* *3v* 
P*(t“T + u * — 


. + vv*—) = p*g y *-^- 


d r ,„dv* 2 du* dv* dw * 

+ dy* ^dy* 3 3^ + 9^* + 3T* )1 
d dv* 3w* d r .du* 9v* 


(3.2b) 


z * component: 

„ ± ,dw* _ *9^* *9^* *9w* x * * dp* 

9f * 3x * dy * dz dz * 

d 9w* 2 du* dv* 9w* 

3 dw* du* d .dv* 9w* 

+ ^ [ ^ + ^ ) 1 + 9 ^ [ ^ + 3 ^ )] 


(3.2c) 
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Energy: 

BT* BT* BT* BT* B BT* 

P * C - b U* h V * h W * ——) = {——— (lc . ——) 

Bt * Bx * By* Bz Bx * 1 Bx * 

Bp* „ Bp* 

+ v *-£— + w *-¥—) + jiO * 

By* Bz 


(3.3a) 


where: 4>* = 2[<f^) ! + (|^)= + 

Bx* By* Bz* Bx* By * 

| |. dw* | 8v* 2 | , 3m* | dw* 2 2 Bu* dv* dw* 2 

dy * Bz* Bz* Bx* 3 Bx* By* Bz* 


(3.3b) 


Note that in the energy equation, p denotes the coefficient of thermal expansion, and d>* 
is called the viscous dissipation function. 

Certain basic assumptions being applied to this problem are: (1) narrow acoustic 
waveguide approximation, yo&«l, (2) no spanwise (z*) dependence and no edge effects, 

^~ = 0, (3) no gravitational field effects, g I *=g y *=g z *=0, (4) no imposed heat 


generation, q" = 0, and (5) fluid is a gas which may be treated as being ideal with 


constant properties (dynamic viscosity, |i, specific heat at constant pressure, c p , and 
thermal conductivity, k c ). The following boundary conditions are also imposed: 


At y* = ±y 0 : u* = 0 

(no slip at wall) 

v* = 0 

(no penetration at wall) 

K,* 

II 

* 

E-h 

(fixed wall temperature, called mean temperature) 

Bu * 

At y* = 0: “ =0 

By* 

(symmetry) 

O 

II 

* 

> 

(“zero-crossover”) 
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The governing equations are transformed into dimensionless equations by scaling 
the dimensional variables using the standard scales for the flow fluctuation quantities in 


an acoustic field. The scaling relations are: 


t = to t*, x = 


(kiln) 


, * y* u* v 

kx*, y = —, u = —, v =- 


(W) 


p - (p*-p.) p - p*-p* t=iizL 

(PJO ’ ' (p m V 0 c)’ AT ref 


where the Mach number is M = —, c is the speed of sound, p m is the mean mass density, 


p m is the mean pressure, and the reference temperature, A T ref , is as yet undefined. 

After implementing the scaling factors, applying the assumptions, and utilizing 
algebraic manipulations in order to simplify the expressions, the terms with order of 
magnitude of ylk 2 « 1 are neglected. The resulting dimensionless equations are: 

Continuity: ^-+ M^-(p«) + ^+ M ^-(pv) + ^ = 0 (3.6) 

at ax ax ay ay 

x - Momentum: 

. aa , AA du. dp ,S 2 .d 2 u (3.7) 

(pM + 1)(— +Mu —+ Mv—) = --^ + (—)—- v > 

at dx dy dx yj dy 

where the boundary layer thickness is 8 = , and the kinematic viscosity is 
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Using order of magnitude analysis and y*ifc 2 « 1, the pressure gradient term is 
the only dominant term in the y-momentum equation (Equation 3.8) which, therefore, 

reduces to = 0. This allows the boundary layer approximation to be made that 
dy 

pressure is dependent only on x and t (p = p(x,t )). 

The remaining governing equations are solved utilizing a regular perturbation 
analysis. The perturbation method assumes a solution of the form: 

/ = /, +e/ 2 +e 2 /j+... (3.10) 
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where £ is a small quantity. In this problem, the appropriate small quantity, the 
perturbation parameter, is the Mach number, M. Therefore, solutions are chosen in the 
forms: 


u = a, + Mu 2 + M 2 u y +... 

4> = +4>i + m $ 2 + M 2 <h + ••• 

M 

p = p, + Mp 2 + A/ 2 p 3 + ... 
v = v, + M\ 2 + Af 2 v 3 + ... 

P = Pr< Mp 2 + M 2 Pi +... 


(3.11) 


Note that the <() solution includes a term of order Af 1 (written O(Af')), since a term of that 
order appears in the energy equation (Equation 3.9). 

These solution forms in Equation 3.11 are substituted back into Equations 3.6, 
3.7, 3.8 and 3.9 which are subsequently expanded algebraically and ordered according to 
the powers of M. The like order terms are grouped, and separate sets of partial 
differential equations are extracted for each order. The most significant solutions are 
those with orders M' 1 and M° (or 0(1)), with the remaining higher order effects not being 
considered in this study. 


Overall, the dimensionless equations to be solved are 


Continuity: 


dp, du. 8v, 

Hr + T± + T' L = 0 

at ox oy 


(3.12) 


x -Momentum: 


3“i _ dPi , 1 

dt dx 2r| 2 By 2 


(3.13) 
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where r| 0 is a gap width parameter based on the viscous penetration depth. 



Energy: 




St 21)1 Pr 3/ p T m 


dt 


(3.14) 


Energy: 0(1): 


3/ 


+ P 


ato,.. 


1 dt 
dp 2 


+ u, 


■ +V 




0 . _ 


dx 1 3 y 2r \l Pr 3y 


1 3 2 <b. . A T ref 

- +2$T m Ec ref 






(3.15) 


The variables, u, <)>, p, v, and p are real functions of time as well as spatial 
variables. In general, a real variable, Z, can be expressed as the real part of a complex 
function as follows: 


Z(x, y, t ) = Re[z(jc, y)e' m ] = ± [ze int + z * *e' int ] (3.16) 

where z(x,y) is the complex amplitude, z**is the complex conjugate of z(x,y ), and 
“n” represents the n" 1 harmonic. Applying this concept to the unknown variables in 
Equations 3.12,3.13, 3.14, and 3.15, it follows that: 
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u x (x,y,t) = R .t[u x {x,y)e i, ] = ^[u x e i ' +u ] **e ~“] 

♦oC^y.O = Re[^ 0 (^,3')e"] = ^10 O <?" +0 O **«'"] 

0 , (x, >,0 = 0 ,0 (*> y) + Re[0 ,2 (*. > , )e 2 " ] = 0 .0 (*. >) +1[0 12 e 2 " + 0 I2 * *e~ 2 " ] 

, (3.17) 

p,(x,>,0 = Re[p I (j:,>)e ,v ] = -[p I e" +p, **<?'"] 

v 1 (*,y,0 = Re[v,(x,>)e"] = ^[v/' +v, **e~ u ] 

JL* 

p l (x,t) = Re[p i (x)e“] = ^[p i e i ' + p, **e~“] 

The variables, 0o, p t , v,, and p y exhibit leading order behavior, whereas 0,, w 2 , 
p 2 , and all others contribute to higher order behavior. These two categories of behavior 
are analyzed in turn in the following sections. The behavior of the pressure distribution, 
however, is particularly significant to this work; therefore, the analysis of p, is treated in a 
separate section following the discussion of the other leading order terms. 

C. LEADING ORDER BEHAVIOR 


The leading order equations to be solved are as follows: 

Continuity: 0(1): —+ ^- = 0 

dt dx By 


x -Momentum: 0(1): ^- = --^- + (—) — — 

dt dx K 2r\l’dy 2 

Energy: O(M'): % = -j- d ? 2 ° + 2»T,Ec(t 


dt 2r)l Pr dy 2 


dt 


(3.18) 

(3.19) 

(3.20) 
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The leading order equations to be solved are now Equations 3.22, 3.23, and 3.24. 
Using the imposed boundary conditions along with the proper equation of state, 
«,, <(> 0 , p P v,, and p } are found. 

The equation of state gives density as a function of pressure and temperature, 
p=p(p,7). A differential element, Ap*, is then given as 
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Ap * =( ^ ViV * +( lr^ ) '- A7 '* <3 - 25) 

Now, 

(^)r* = pr, A p* = p*-p m = pp m V 0 c (3.26) 

Also, 

A!T* = r * -r m = §AT ref and Ap* = p *-p m = pp m M (3.27) 


From the definition of the coefficient of thermal expansion, B = —— (^—) ,, it follows 

p* BT* p 


(|^V = "Pp* = -Pp m d + Mp) (3.28) 

Utilizing all of these relations. Equation 3.25 becomes 

PP„A* = ^Pp m V 0 c- p mPm (l + Mp^AT ref (3.29) 

Substituting in the perturbation solution forms from Equation 3.11, the following 
equation results for pi: 

P, ='VP,-P m 7' m (y-lX>o (3.30) 

Using the complex notation established in Equation 3.17, it results that: 

Pi (3.31) 

For ease of following this discussion, the final leading order equations to be solved and 
their corresponding boundary conditions are summarized below: 

„ . . 3v. _ Bu. 

Continuity: — L = -,p ] _ _ L (3.32) 

ay dx 
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(3.33) 

(3.34) 

(3.35) 


The density, p,, comes directly from the solution of the equation of state 
(Equation 3.35). That leaves three equations and four unknowns which, at first glance, 
are impossible to solve. However, a fourth equation is obtained when applying one of the 
boundary conditions to the continuity equation (Equation 3.32). When applying the 
condition that v, =0 at y=0, an equation for p,(x) results. Then, using the appropriate 
boundary conditions for p l , the function, /?,(*), is determined explicitly. The expression 
for p t (x) is then utilized in the remaining equations to solve explicitly for the remaining 
unknowns, «,,<!>„, p,, andv, . 

Using standard methods of solving partial differential equations (such as 
separation of variables). Equations 3.32, 3.33, 3.34, and 3.35 were solved. Some 
simplifying substitutions were also made. Since the reference temperature is defined to 
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be A T ref = , the Eckert number, as defined in Equation 3.9, is simply unity. Also, for 

2c p 

an ideal gas, the quantity, (3 m 7 m is approximately unity. Thus, the leading order solutions 
are given by the following expressions: 


dx cosh((l + i)n 0 ) 


(3.36) 


t- , cosh((l + i)ri 0 y->/Pr) 

Q> 0 (x,y) = 2 Pl (x)[l -— r—A 

cosh((l + i)q 0 yvPr) 


(3.37) 


cosh((l + t)r| 0 vPr) 


(3.38) 


v, (x,y) = ip { (x)[(l - y) + (y - 1)( 


tanh(( 1 + i)T| 0 VPr) 

(i+/)r| 0 Vp7 


sinh((l + /)q 0 y-\/Pr) . J 2 p, 

(l + i)T|o'v/^cosh((l+/)r| 0 VPr) dx 2 ^ 

+ ( tanh((l + Qq 0 ) _ sinh((l+ Qq 0 y) 

(l + /)ri 0 (l + Ori 0 cosh((l + i)ri 0 ) 


(3.39) 


D. THE PRESSURE DISTRIBUTION 


Applying the boundary condition that v, = 0 at y=0 to the continuity equation 
(Equation 3.32), an ordinary second order differential equation for p x (x) results as 
follows: 


d 1 p. 9— 

—r L + w p.(x) = 0 where 

dx 

m 2 __ 0+ *)n 0 vPr 

^ tanh((l + t)n 0 ) 

d + OHo 


(3.40) 
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The boundary conditions necessary for determining p i (x) are found in the geometry of 
the problem. As shown in Figure 3.1, the plates are flush with the source of the acoustic 
wave at x=0. This means that the pressure amplitude is at a peak (magnitude of one), and 
the velocity amplitude, and therefore the pressure gradient, is at a node (magnitude of 
zero). Thus, the appropriate boundary conditions needed for determining the pressure 

distribution are at x=0, p, = i and —= 0. Solving for p, (x) and applying these 

ax 

boundary conditions, the result is 

p,(x) = i'cos(mx) (3.41) 

This expression is then substituted back into Equations 3.36, 3.37, 3.38, and 3.39 to give 
explicit forms of w,, <J> 0 ,p,, and v, respectively. 

The two cases studied in this research are the wide gap and arbitrary gap cases 
which differ bases on the width of separation between the plates as compared to the 
Stokes’s boundary layer thickness as discussed in Chapter I. The basic assumption that is 
made in the wide gap case is that, in the pressure distribution of Equation 3.41, the 
parameter “m” is assumed to be unity since m~\. This simplifies the problem 
significantly, since the pressure distribution just varies as cos(x) and is independent of any 
other parameter. For the arbitrary gap case, “m”, is left as a function of T| 0 , y, and Pr as 
shown in Equation 3.40, and the pressure distribution becomes dependent upon those 
three parameters as well as on the variable, x. 
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E. HIGHER ORDER BEHAVIOR 


The higher order velocity, density, and pressure terms, w 2 , p 2 , v 2 and p 2 , and 
higher will not be considered in this study; however, the higher order temperature term, 
<J>i, is significant and will be considered in this section. The equation to be solved in this 
consideration is the 0(1) energy equation: 


Energy: 0(1): 


fro ! u fro , v fr _ 1 ^ 2 4>i 

dt 1 dx 1 dy 2r\l Pr fr 2 


+2pr m Ec— 




(3.42) 


j 

Again applying the equality, ~~~ = M 2 (— -), and subsequently neglecting terms of 

T* 2 

nt 

order A/ 2 , Equation 3.42 further reduces to 


Energy: 0(1): 


frl + D fr<L + M fro , v fro = 1 dfr. , 

dt 1 dt 1 dx 1 dy 2r\l Pr dy 2 


dt dx fr 


(3-43) 


Equation 3.43 is not complete; it is missing one term. The extra term comes from 
expanding the coefficient of thermal expansion, p, into its perturbation form. By 
definition, 


Q = _ 1 r fr*^ __ 1 r/fr*. . AT * , 

^ p* 97* * p* 5r* m 3r* 2 m 


(3.44) 


where: 


A T* = T*-T m =<>A T re 



Utilizing the appropriate scaling factors for p* and T* from Equation 3.5, and the 
corresponding perturbation solution forms from Equation 3.11, the expanded equation for 


pT* becomes 


Pr* = P . {1 - M[p, + -1 ] + 0(M ! )+... 1 

p*m rmr/n 


(3.45) 


Since 


T d 2 p * 

m / r \ _ 


p.p. w 


( , ) = 2 for an ideal gas, and 


c 2 y -1 
2c r " “T~ 

p m 


, Equation 3.45 simplifies 


to 


Pr* = P.r.(i-M[p, +IL!-)*,+0(M J )+...) 


(3.46) 


Substituting the expanded P7 1 * expression into the energy perturbation equation, an extra 

0(1) term comes out of the 2£cpr *(—)— term from the O(M) portion of P T 

M dt 

multiplying with (1/M). Therefore, the extra 0(1) energy term is 

~2P m f m [pi + . making the complete 0(1) energy equation as follows: 

dt 2 at 


Energy: 0(1): 

^L +f>i ^ + Ui ^± +yi ^ = _^_^± + 2KTE A 

dt dt dx dy 2r| 2 Pr dy 2 dt 

+ 2P A to h 4L + S(*!L)._ 2 p. r j Pi ai + (lzi)*.3EL ] 


(3.47) 


dx r) 2 0 dy 


dt 


dt 


Substituting the relations from Equation 3.17 into Equation 3.47, yields the 
following complex partial differential equation: 
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Energy: 0(1): 


,2 it , rz a „2i< , — a<t>0 2it . — 9*0 _2/r 


Re[2i<|> 12 e + ip,<|) 0 e 2 " + m 


i a '« + v . 

ox 


dy 


■.*»!_ 1 3 >10 

2rio Pr dy 2 


+ Re[ 


1 - *12 „2it , ,, r „R t „2,'t , 3 m, n2 _2i( 


2^1 o Pr ay 2 

- dPx 


-e +4Ec$ m T m ip 2 e ' + —j(-jp-) 2 e 

T| 0 ay 


(3.48) 


+2£cp m r m i7 1 -^-* 2 " - 2p m r m /p, p., 2 " - p m r m /(Y - i)* 0 p,* 2 " 1 


Two equations may be extracted from Equation 3.48. One yields information about the 
first harmonic of <J),, which is the real variable, <|> 10 ; solution of the other equation yields 

the second harmonic, <|> 12 , which is complex. The determination of 0 12 is not of interest 
in this study since it does not contain a steady temperature component; however, (|) 10 is 

the time-averaged temperature. Since it is of great importance to this research, the time- 
averaged temperature, <() 10 , is considered in detail in the next section. 


F. THE TIME-AVERAGED TEMPERATURE, <t> 10 


In Equation 3.48, there is one term which obviously has no time dependence. 


( 


i a 2 <t> 


10 V 


2ilo Pr 3y 2 


); the remaining terms are multiplied by e 2 ". This supposed time 


dependence would imply a simple solution for <J> 10 since it would mean that 


d 2 <t>,o 

ay 2 


= 0 . 


However, upon close inspection of the e 2 " terms, it is noted that most of them have a 
steady (time independent) term embedded in them. Taking the time average of Equation 
3.48, the only terms which drop out (average to zero) are 2/<|> 12 e 2 ", 
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however, the calculus and algebraic manipulations involved in solving the equation by 
hand are extremely cumbersome. For example, evaluating the time average of the first 
term, < ip,<|) 0 >, is accomplished using the definition of time average from Equation 3.50 


and substituting into it the expressions for p, and <|> 0 from Equations 3.37, 3.38, and 
3.41 as follows: 
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<P,i*. >= \ Re[p,H’)»o **] = \ Re(tcos(ffu:)[l + COSh((l + ^ ♦ 

2 2 cosh((l + i)n 0 vPr) 


(-i)[2i‘cos(mx)(l - 


cosh((l + /)n 0 yVPr) 
cosh((l + i')n 0 VPr) 


•)]**} 


where m = 


(1+ xiaj^n^F) 
_(1 + OnoVPr 


tl ( tanh((l + i)n 0 ) 
0 + i)Ho 


(3.51) 


Once the complex conjugate of <J> 0 is correctly computed, the real part of the entire 
expression must be evaluated without error. This establishes the first time average term, 
but there are six more which must be obtained in a similar fashion. The summation of all 
seven resulting time average terms must then be integrated twice, applying appropriate 
boundary conditions, to obtain an expression for <f) 10 . The vast amount of careful 
mathematical manipulation and computation required to solve for <j) l0 points directly to 

the use of a powerful symbolic mathematical software tool. Thus, MAPLE, just such a 
computer software package, is utilized to solve for <|> 10 , given all of the other expressions 

for m, , <|> 0 , p,, v,, and p ] . <|> 10 is found first for the simpler wide gap case and then for 
the arbitrary gap case. The resulting expressions for <(> 10 are much too lengthy to include 
explicitly in this paper; however, the temperature distributions described by <|> l0 are 
plotted and discussed in Chapter IV. 
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IV. RESULTS AND DISCUSSION 


The results of this work are comprised of the solutions to the governing equations 
as given in Chapter in along with plots of these solutions. These plots are presented and 


discussed in the following 


thermal penetration depth is 


sections for different values of n 0 VPr = — where the 

5„ 


2 a ^ 

8 a =, — = -p=. The parameter, Ti 0 VPr, has been found 
V to vPr 


to be the more appropriate gap width parameter than T | 0 in dealing with thermal effects in 
such oscillatory flow problems. Since the form of the leading order pressure distribution. 
Pi (*), is crucial in determining the behavior of the remaining leading order solutions and 
the higher order solutions, the pressure distribution results are analyzed first. Then, the 
leading order velocity, temperature, and density solutions and the higher carder time- 
averaged temperature term are presented. Finally, the recovery factor, suitably defined 
for this model problem, is presented, plotted and discussed. 


A. THE PRESSURE DISTRIBUTION 


The pressure distribution, /?,(*), for this model problem is found, as given in 
Chapter III, as 


Pi CM) = Re[p, (*)£"], where p t (x) = i cos(mx) and 


m = 


D+( ,-ixH!«a±a^) M 

_(1 ~t~ f)f|Q>/Pr 


0 ^ tanh((l + i')n 0 ) 

O + OHo 


(4.1) 
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Note that the eigenvalue “m” is in general complex; however, with increasing gap width, 
the parameter “m” approaches a value of one. Therefore, in the idealized wide gap case, 
it is assumed that m= 1, and the pressure distribution is a simple cosine function, i.e. 
p x (x) = /cosx and p x (x,t ) = -cosxsin/. In general, for the arbitrary gap case, “m” is a 

complex number, m = a - bi. As gap width (or ti 0 ) increases, “a” reduces to a value of 
one and “b” to a value of zero, approaching the results of the idealized wide gap case. 
For all of the trials calculated, the values of “m” are given in Table 4.1. The value of “m” 
also depends on the gas properties, y and Pr. Thus, for the same T | 0 value or gap width, 
the magnitude of “m” will vary for different gases. The ideal gas utilized in this work is 
air; therefore, the property values used are y=l .4 and Pr=0.7. 


Table 4 

kl Values of “m” in Pressure Distribution. 

o 

m 

M 

1 

1.3047-0.5939/ 

1.443 

1.5 

1.1814-0.3165/ 

1.223 

2 

1.1405-0.2055/ 

1.159 

3 

1.1002-0.1229/ 

1.107 

4 

1.0762-0.0884/ 

1.080 

5 

1.0613-0.0688/ 

1.064 


Two plots of the pressure distribution, p x (x,t), are given in Figure 4.1, for T| 0 VPr 
values of 1.5 and 5. Considering the ri 0 VPr =5 case, the pressure distribution closely 
resembles the idealized wide gap case of p x (x,t) = - cos*sin /. As gap width decreases, 
the amplitude of the cosine function increases, and the waveform is considerably distorted 
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as illustrated by the T| 0 VPr =1.5 plot. This variation in amplitude is clearly observed in 

TC 

Figure 4.2 where p,(x,t) is plotted at two particular times t=— and t=n for r^-v/Pr values 

of 1.5, 2, 3 and 5. To emphasize the distortion as compared to the idealized wide gap 
case, the function, p { (x,t) = -cosxsinf is also plotted in Figure 4.2. Based on the 
calculated values of “m” and the plots of the pressure distributions, the gap width 
corresponding to r| 0 VPr =5 is determined to the sufficiently close to the idealized wide 
gap case to be considered as a “wide gap”. Again as r\ 0 decreases into the “narrow gap” 
region (r| 0 VPr <5), the pressure distribution significantly deviates from the simple cos(.x) 
behavior of the idealized wide gap case. 

B. LEADING ORDER BEHAVIOR 


The leading order first harmonic solutions for the primary variables (velocity 
components, temperature and density) are given as follows (as discussed in Chapter HI): 
M,(x,y,r) = Re[n,(*,yy'], ty 0 (x,y,t) = Re[<j> 0 (;t,;y)e"], p,(jt,y,r) = Re[p,(jt,y)e"], and 


Vi(x,y,0 = Re[v, (*,)>>"], where 


u,(x,y) = - c °sh((l + »),>’) 

dx cosh((l + t)r| 0 ) 


<t> 0 (^,V) = 2^» 1 (at)[1 — 


cosh((l + /)T| 0 yVPr) 
cosh((l + i)ri 0 VPr) J 


Pi(*>;y) = Pi(*)[Y — (T — 0(1 — 


cosh((l + i)n 0 yVPr) 
cosh((l + /)r| 0 VPr) ^ 


(4.2) 

(4.3) 


(4.4) 


32 





v, (x, y ) = i Pi (*)[( 1 - y) + (Y - 1)( 


tanh(( 1 + pT) 0 VPr) 

(l + Otlo^Pr 


- sinh((l + ,XWft) j^L [( 0 (4 . 5) 

(l+i)T| 0 Vfr cosh((l+T)ii 0 VPr) dx 2 

+ ( tanh((l + Qri 0 ) _ sinh((l+ Qq 0 y) 

(l + /)Tlo (l + On 0 cosh((l + i)r| 0 ) J 

These quantities, as opposed to the pressure distribution, depend on both spatial 
variables, x and y. Their x dependence is given by either cos(mx) or sin(mx) which stems 
from their dependence on the pressure distribution or the pressure gradient. Because of 
this periodicity in x, all of the information about the x dependence of a particular leading 
order distribution is contained in an interval of one half of the wavelength (i.e. 
0<*<7t). In the y direction, the quantities u x , <J> 0 , and p, are symmetric about the 
centerline of the channel (y= 0 ), whereas the velocity in the y direction, v,, is 
antisymmetric about the centerline. Therefore, the results for all of these distributions are 
determined and presented for one half of the channel only, specifically for 0 < y < 1 . The 

plots of the leading order density distribution, p, (x,y,t), are not included in the discussion 
since the behavior of the density function is not of direct interest to this work. 

The velocity distribution, u x (x,y,t), is given in Figure 4.3 for x=n and q 0 VPr 

values of 1.5 and 5. u x (x,y,t) depends on the pressure gradient; therefore, it varies in x as 
sin(mx). This distribution, as observed in Equation 4 . 2 , varies in y as 

Re[l-- 7 -— '^°^ 3. Figure 4.4 shows u x (x,y,t) for x=n at two instants in time, t=— 

cosh((l + i)n 0 ) 2 

and t=n for r | 0 VPr = 1.5, 2 , 3, and 5. The distributions clearly show that the velocity 


vanishes at the wall (y=l) and thereby satisfies the no slip condition. Also, because 
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Ui(x,y,t) is symmetric about the centerline of the channel, the gradient there is equal to 
zero as shown by the zero slope at y=0. For wide gaps or large values of T| 0 , the velocity 
component, u u increases in magnitude from zero at the wall (y=l) to a maximum value 
and then decreases slightly to level off with a slope of zero at the centerline (y=0). As the 
gap narrows or "Ho decreases, the w, velocity value increases, and the location of the peak 
velocity in the channel shifts toward the centerline. In fact, for T| 0 VPr =1.5, the 

maximum velocity value occurs at the centerline, and the distribution is more parabolic in 
shape (similar to that of Poiseuille flow). These trends are clearly observed in the curves 

7t 

for both snapshots in time, t=— and t=n, in Figure 4.4. 

The leading order temperature solution, § Q (x,y,t), is illustrated in Figure 4.5 for 
x=n and T| 0 VPr = 1.5 and 5. The x dependence is a cosine function, icos(mx), which is 
caused by the direct relationship of <f> 0 to the pressure distribution. The y dependence of 
, . e n r .„ cosh((l + On ft yVPrT, 

<t>o is of the form Re[i(l--— - ,— )] as observed in Equation 4.3. In Figure 

cosh((l + i)r| 0 VPr) 6 

4.6, where <J> 0 is plotted for x=n, at two instants in time, t=— and t=K, and for T) 0 VPr 

values of 1.5, 2, 3, and 5, the dependence of <J> 0 on y and T| 0 is displayed, and its similarity 
to that of the u t velocity distribution is noted. The boundary conditions of no slip at the 
wall and symmetry at the centerline are clearly met as <t> 0 =0 at y=l and its slope vanishes 
at y= 0. Considering the curves for time t=n, for large gap widths or T | 0 values, the 
temperature, <t> 0 , increases in magnitude from zero at the wall to a peak and then decreases 
to the centerline temperature value which is approached at zero slope. Again, as the gap 
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width or r| 0 value decreases, the peak shifts toward the centerline and its magnitude 
increases. As with the m, velocity component, the shape of this temperature distribution 
becomes more parabolic as the gap narrows or T | 0 decreases. The dependence of the 

K 

temperature, <t> 0 , on gap width or r| 0 is not accurately portrayed in the plots for time, t=—. 

This is explained by the nature of the time dependence of the temperature distribution, 
<t>o(*,y,0- A phase lag occurs in time, and the amount of this phase lag depends on the 

K 

value of r| 0 ; therefore, the temperature distributions at the instant of time, t =—, for the 

different gap widths do not correspond to the same time positions in their respective 
periods. 

The velocity component in the y direction, v x (x,y,t), is presented in Figure 4.7 for 
x=n and r| 0 VPr = 1.5 and 5. From Equation 4.5, it is observed that Vi(x,y,t) depends on 

both the pressure distribution and its second derivative. Therefore, the x dependence is a 
cosine function, cos(mx). The y dependence of v,(x,y,t) is a combination of linear terms, 
hyperbolic sine terms and constants, resulting in a function which is antisymmetric about 


K 

y= 0. In Figure 4.8, v,(x,y,f) is plotted for x=K, at two instants in time, t=— and t=n, and 

2 


for T| 0 VPr values of 1.5, 2, 3, and 5. The boundary conditions of no penetration at the 
wall and “zero crossover” at the centerline for an antisymmetric distribution are satisfied 
as v,=0 at y=\ and y=0 respectively. The behavior change in v, as gap width or r| 0 
changes is clearly illustrated by the *=tc curves. For large gap widths, the y component of 
the velocity, v,, increases in magnitude from zero at the centerline to a peak near the wall 
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and then decreases rapidly in magnitude back to zero at the wall. An inflection point is 
observed between the centerline and the peak value for wide gaps. As the gap narrows or 
Ho decreases, the peak increases in magnitude and shifts away from the wall. Also, the 
shape of the distribution becomes more parabolic for narrower gaps (i.e. the inflection 
point observed in the wide gap trials disappears as Ho decreases). The trend cannot be 

K 

distinguished in the t=— plots for reasons relating to the nature of the time dependence of 

£ 

v,(jc,y,0 which are similar to those discussed in the <t> 0 (*,;y,0 results. 

C. THE TIME-AVERAGED TEMPERATURE 


The time-averaged temperature distribution, (|)|o(jc,y), is shown in Figure 4.9 for 
H 0 VPr = 1.5 and 5. From the partial differential equation for <t>, 0 as presented in Chapter 
El and again here for convenience, 


< ip,<t> 0 > + < u, > + < v, ^2- > 


1 9 2 <|) 


,0. Jc d«, x2 


Bx 

— dP\ 


+ <■ 


By 2r\l Pr By 2 Ho 3y 


■(rzr^r > 


(4.6) 


+ <2Ec$ m T m u l -£>-<2$ m T m ip l p l >-<P m r ra /(Y-l)<() 0 p, > 


the x dependence is determined by the various time-averaged terms in the forcing 
function. The terms involve the multiplication of two of the leading order first harmonic 
solutions of the primary variables, <J» 0 , p, , and v,, or their derivatives. Substituting the 
appropriate distributions from Equations 4.2, 4.3, 4.4, and 4.5 into Equation 4.6, the x 
dependence of <|>jo is observed to stem from terms where the pressure distribution is 
multiplied by itself as well as terms where the pressure gradient is multiplied by itself. 
Therefore, the resulting time-averaged temperature distribution varies in x as a 
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combination of sin 2 (mx) and cos \mx) terms. In Figure 4.10, the time-averaged 

temperature, <M*oO, is shown for x=0, x= — and x=n and for T| 0 VPr values of 1.5, 2, 3, 

2 

and 5. The temperature, <t»i 0 , vanishes at the wall (y=l), satisfying the no slip boundary 
condition; and its gradient,= 0 at the centerline (y=0), thereby satisfying the 

symmetry condition. From the x=0 curves, the variation in <j>| 0 based on T| 0 is clearly 
observed. For wide gaps or large rj 0 values, the time-averaged temperature increases from 
zero at the wall to a maximum value fairly close to the wall and then decreases gradually, 
leveling off to a zero slope at the centerline. Thus, for wide gaps most of the temperature 
variation occurs near the wall. As the gap size becomes smaller or T| 0 decreases, the 
magnitude of the time-averaged temperature gradient decreases, and the temperature 
variation is spread more and more evenly across the channel width. Further discussion of 
the affect on the time-averaged temperature gradient of varying the gap width is presented 
in the next section on the recovery factor. The time-averaged temperature distribution, 
is not symmetric in x. This is illustrated by the fact that the plots for x=K do not 
overlap with the curves for x=0. 

D. THE RECOVERY FACTOR 

More important than the results for the time-averaged temperature distribution, 
<|>io, presented above, it is the gradient of this temperature at the wall which allows key 
deductions to be made about the time-averaged heat transfer effects in this problem. In 
order to interpret these heat transfer results in a more conventional manner, they are 
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presented in terms of a suitably defined recovery factor. The recovery factor, r, is the 
ratio of recovered thermal energy at an adiabatic surface to the fluid’s kinetic energy as 
discussed in Chapter II. The form of “r” utilized in this research conforms to convention 
and is given by 

v o 


T m is the fixed temperature of the upper plate. T uw is the adiabatic wall temperature and it 
is equal to the value of T* at y*=-y 0 when the wall there is adiabatic. The adiabatic 

8r* 

condition is that the temperature gradient, -— = 0 at the wall in question, y*=-y 0 . In the 

By * 

present problem, T aw would correspond to the time-averaged temperature achieved by a 
wall if it were made adiabatic. Let the corresponding dimensionless temperature for this 
adiabatic wall be denoted in general as \|/. Then it is the corresponding time-averaged 
temperature, V|/ 10 , at the adiabatic wall that is of interest here. Expressing the recovery 

factor and the adiabatic conditions in terms of the dimensionless parameters yields: 


-jJLAT ref =0 at y = -1, 

r aw =\ l / 10 Ar re/ +r m at y = -i, 
c p (M>io AT nf +T m — T m ) \|f w AT r 


r-=¥. o> 


or r =\\f 10 (y = —1) 

where y ]0 is the dimensionless time-averaged temperature obtained in the adiabatic wall 


47 



In order to find “r” in this manner, however, it is necessary to additionally solve 
for the temperature distribution, \jt l0 , in the adiabatic wall case. To avoid this extra time- 

consuming step, a relationship will be proven between the value of the time-averaged 
temperature, \|t 10 , at an adiabatic wall and the value of the time-averaged temperature 


gradient, , at the wall in the fixed temperature case that is the subject of this 


research. The two distributions appear to have approximate shapes as given in Figure 
4.11. 



Figure 4.11 Approximate Temperature Distributions for Two Cases: 
(a) Fixed Wall Temperature and (b) Adiabatic Wall Temperature 
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The recovery factor has already been shown in Equation 4.8 to be equal to 

Vio(:y = “!) • 11 will now be proven that \j/ 10 (y = -1) = 2^-(y = -1). From the nature 

dy 

of the partial differential equation for (J) l0 in Equation 4.6, it is known that the y 
dependence of <f> 10 can be expressed as 

<M>0 =/OO + C,y + C 2 (4.9) 

where from the nature of the forcing function in Equation 4.6, it can be shown that/(y) is 
an even function. Now, applying the appropriate boundary conditions at the walls for the 
present problem yields 

<>,o(D = 0 = /(l) + C,+C 2 and 

<t>io( - l) = 0 = /(—1)-C, + C 2 (4 - 10) 

Since/(y) is an even function, /(l)=/(-l), and the constants of integration are determined 

to be C, = 0 and C 2 = -/(1). Thus, the resulting expressions for <|> 10 and -ft 10 are 

dy 

<l>io (y) = f(y)-fV) (4.ii) 

(4.i2) 

Similarly, for the adiabatic wall problem, the time-averaged temperature has the 
same governing equation since the primary leading order variables are the same. Then, 
\|r 10 is expressed as the same function of y but with different constants of integration as 
follows: 

Vu>(y) = f(y)+D iy +D 2 (4.13) 
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Applying the appropriate boundary conditions for this case, which are \j/ lo =0 at 


y=l and 10 =0 aty=-l, yields 
By 


Vio(l) = ° = /(!) + £>, +D 2 and 
^(1) = 0 = /'(1) + D. 


(4.14) 


Thus, the constants are found to be D, = -/'(-l) and D 2 = and the 

resulting function is found to be 


v 10 0» = f(y) - f'i-Dy + /'(-I) - /(l) 


(4.15) 


Evaluating the adiabatic wall temperature, y 10 at y=-l, yields 


V ,o(-D = /(-l) + /'(-l) + /'(-l) - /(l) 


(4.16) 


Again, noting that/(y) is an even function and using the properties of even functions, the 
previous equation is reduced to 


\|/ 10 (-l) = 2/'(-l) 


(4.17) 


Also, from Equation 4.12, it is seen that -^-(y = -1) = /'(-1). Therefore, 

By 


\K, 0 (-l) = r = 2^(y = -l) 

By 


(4.18) 


This proof allows “r” to be defined and calculated in terms of the temperature 
distribution, <|> 10 , which has been determined in the present problem without having to 


redo all of the work necessary to solve the adiabatic wall case. The gradient, —, is 

By 


easily obtained and evaluated at the wall and then doubled. Since f(y) is an even function, 
/'(y)is an odd function; therefore, it is known that/'(-l) = -/'(!). The domain of 




analysis is 0 < y < 1; therefore, for convenience, the recovery factor calculated in this 
research is determined by 

r = _2 ^r L 0’ = 1) (4.19) 

dy 

The resulting values of “r” may be positive or negative. Positive values of “r” 
correspond with the physical situation of the wall being heated by the fluid. This 
phenomenon of the walls being heated by viscous dissipation is observed in high speed 
mean flows such as in the example cases presented in Chapter II. Negative values of the 
recovery factor indicate that the direction of heat flow is away from the wall as shown in 
Equation 4.19. In other words, the wall is actually being cooled by the fluid in those 
instances. This cooling potential is a very interesting attribute of high speed oscillatory 
flows, and it is this phenomenon which is the underlying principle that makes 
thermoacoustic refrigeration possible. Thermoacoustic heating and cooling effects caused 
by thermoacoustic streaming in a resonance tube were predicted theoretically and then 
validated experimentally by Merkli and Thomann (1975). A subsequent review of the 
topic was provided by Rott (1980). Experimental results were also obtained by Wheatley, 
et al. (1983) which verified the existence of this cooling effect in a stack of plates placed 
in an acoustically resonant tube. 

The recovery factor, r, is calculated using Equation 4.19 and then plotted for each 
trial run for this problem. The results are illustrated in Figure 4.12 for one half of a 
wavelength, 0 < x < n , and for T) 0 VPr =1.5, 2, 3, and 5. For each gap width examined, 
the recovery factor takes on both positive and negative values which corresponds to 
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heating and cooling of the channel wall respectively. In general, at jc= 0, the recovery 
factor is a positive local maximum. As x increases, “r” decreases for each gap width until 
it reaches a minimum value which is negative. The recovery factor then increases as x 
increases until it reaches its maximum positive value before again beginning to decrease. 

As observed in Figure 4.12, the recovery factor is generally higher in magnitude 
for all x values plotted and the curve is steeper (i.e. the gradients are higher in magnitude) 
for a wide gap than for a narrow gap. The fact that the “r” values are higher for wide 
gaps than for narrow gaps corresponds with the higher time-averaged temperature 
gradients observed near the wall (y=l) for wide gaps in Figure 4.10. According to 
Fourier’s law as well as the meaning of the recovery factor as described earlier, these 
observations would imply greater heat fluxes for greater gap widths. This apparent 
implication does not agree, however, with what has been experimentally observed. 
Therefore, in order to better understand physically what “r” is measuring in this model 
problem of oscillatory flow, the equation for the recovery factor utilized in this work. 
Equation 4.19, is expanded by substituting back into it its dimensional quantities. From 
Equation 4.19, the recovery factor is proportional to the dimensionless time-averaged 
temperature gradient as follows: 


r oc 


*►■0 

dy 


(4.20) 


Using the appropriate scaling relations involving <)) l0 and y from Equations 3.4 and 3.5, 
the relationship in Equation 4.20 becomes 
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(4.21) 


r oc ->a .j£Q * 

X v 3y* 

Manipulating the expression in Equation 4.21 in order to relate the recovery factor to the 
physical parameters which are more appropriate for this oscillatory flow problem 

produces 


4xf) ** 

c 2c„ 


(4.22) 


where =Ti 0 VPr is again the gap width parameter utilized in this work. Substituting 
6 a 

y 

the relations, M = — and c 2 =yRT m , into Equation 4.22, the expression for the 
c 

recovery factor becomes 


( Sq } ^C%) = 'Ho^Pr d( r ' n X.) 

*(%) ~ M 2 ( Y^ 3(>X) 
2 2 


(4.23) 


From Equation 4.23, the apparent contradiction to physics observed in Figure 4.12 
is explained. It is clear from Equation 4.23 that for increasing values of t| 0 VPr, an 

increase in the magnitude of “r” does not automatically imply an increase in the time- 
averaged temperature gradient on the right hand side (RHS). In order to distinguish 
clearly the effect on the time-averaged temperature gradient as the gap width (or T| 0 VPr) 
changes, the recovery factor, r, must be scaled by T3 0 VPr as follows: 
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r 


(4.24) 


_ 1 3( r 'vO 

n,VF “ ^izl) S(’X) 

2 

Therefore, for a given Mach number, M, and for a particular gas (i.e. y fixed), the 
r 

quantity,-, in Equation 4.24, carries the physical information about how the time- 

r| 0 VPr 

averaged temperature gradient (and thus the heat flux) varies with different gap widths (or 
values of ri 0 VPr). 

The recovery factor, scaled by T| 0 VPr , is plotted in Figure 4.13 for r| 0 VPr = 1.5, 
2, 3, and 5 and illustrates the true heating and cooling effects at the wall (y=l). For a 

wide gap as compared to a narrow gap, the “scaled recovery factor” (— r - T =) is lower in 

■HoVPr 

magnitude when x=0 and remains lower until approximately x= 0.271. At this point all of 
the curves cross so that the scaled recovery factor is slightly higher in magnitude for wide 

gaps until it crosses zero. During the initial cooling portion of the curve (—^=<0), the 

tloVPr 

wide gap magnitudes are lower again until the curves turn at their respective minima. 

Once the heating region is reached again (—4=>0) for all of the curves, the scaled 

iloVPr 

recovery factor is again lower in magnitude for wide gaps. The —curves are flatter 

tloVPr 

(i.e. gradients are lower in magnitude) for wide gaps for generally all values of- j=. 

Ti 0 VPr 

Considering the portion of the curves where the scaled recovery factor is negative, the 


55 





following observations are noted. For r| 0 VPr=5, the zero crossings of the scaled 
recovery factor occur at approximately x=0.36tt and Jt=0.587t. Thus, the length of the x- 

V 

interval for which — —j= is negative is approximately 0.22k. As the gap width narrows 

or T| 0 decreases, the average magnitude of the scaled recovery factor decreases and the 
magnitude of its gradient increases. This results in a shift of both zero crossings to lower 
x values as T| 0 decreases as well as a shortening of the ^-interval over which the scaled 
recovery factor is negative. For example, when r| 0 VPr =1.5, the zero crossings of the 

scaled recovery factor occur at approximately x=0.33n and x= 0 . 467 t, resulting in an x- 
r 

interval where the- j= value is negative of approximately 0.1371. 

TloVPr 

The higher magnitudes of the positive scaled recovery factor for narrower gaps 
translate into higher heat fluxes at the wall of the channel for those x locations. Just as in 
the previous discussion concerning the conventional recovery factor, r, for the x-intervals 
where the scaled recovery factor is positive, the channel wall is being heated by the fluid 

through viscous dissipation. For the ^-interval where —— f — is negative, the channel 

TlovPr 

wall is actually being cooled by the oscillating fluid. The gradient of —is also 

Ti 0 VPr 

important because it determines the heat pumping capacity available. More specifically, 

for the same magnitude of-(or cooling effect), a higher gradient equates to more 

TloVPr 

heat pumping capacity available for the same length of plates. 
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This cooling effect, produced in zero-mean oscillatory flows, is utilized in 
thermoacoustic refrigeration applications. The stack of the thermoacoustic refrigerator is 
positioned in a tube where part of the stack is in a region of x values where the scaled 
recovery factor is positive, and the remaining portion of the stack is in a region where 

T f 

- 7 ~ is negative. The end where cooling occurs (negative -p= values) is where 

r| 0 vPr 

the refrigerating effect is produced. The heat generated at the other end of the stack 
(where the scaled recovery factor is positive) is removed via a heat sink. The stack region 
of a typical thermoacoustic refrigerator is very short compared to the wavelength of the 
acoustic wave producing the oscillatory flow; therefore, an analysis of the scaled recovery 


factor in the region of x where 


r 

TIoa/P^ 


just becomes negative is important for this 


application. 


Comparing both the magnitudes of- j= 

T] 0 VPr 


and the gradients in the region 


of x from approximately 0.33 tc to approximately 0.4rc in Figure 4.13, it is observed that 
the highest magnitude and gradient of the scaled recovery factor occur for the case where 
T| 0 VPr =2. Therefore, for a stack place with its cooling end in that x region, the heat flux 


away from the wall and the heat pumping capacity are both greater for a gap width 
corresponding to T| 0 VPr =2 than for any of the other gap widths examined. 
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V. CONCLUSION 


This work involved finding the analytical solution to a fundamental problem in 
order to study the time-averaged heat transfer effects induced by high speed zero-mean 
internal oscillatory flows. The governing equations of continuity, momentum, energy, 
and state were established for the model problem which consisted of an ideal-gas host 
fluid between two parallel plates excited into a zero-mean oscillatory flow by a high- 
intensity standing acoustic wave. 

Introducing the proper scaling relations, the governing equations were 
transformed into dimensionless equations; then, order of magnitude analyses were 
conducted to simplify the resulting equations. Employing perturbation solution methods 
with the Mach number, M, used as a perturbation parameter, the equations were expanded 
and separated based on the orders of M. The most significant order equations were then 
analytically solved using complex variable methods for the unknown velocity, density, 
temperature, and pressure terms. The final aim was to determine the time-averaged 
temperature distribution for various gap widths in order to analyze the heat transfer 
effects induced by the zero-mean oscillatory flow. 

The analysis was completed using air as the host fluid for gap widths 
corresponding to various T| 0 VPr values up to T| 0 VPr =5, which was determined to be the 
lower limit for consideration as a “wide gap”. The leading order primary variables were 
presented graphically and discussed in detail in Chapter IV along with the time-averaged 
temperature and corresponding recovery factors. According to the interpretation of these 


59 


results, using air as the host fluid, the most suitable choice of gap width in terms of the 
cooling effects applicable to thermoacoustic refrigeration corresponded to a value of 
r| 0 VPr=2. 

This was a fundamental study, hoped to be a foundation upon which to be built, in 
the area of recovery factors in zero-mean internal oscillatory flows. Further work is 
suggested such as determining the most suitable gas to be used as a host fluid (i.e. to 
maximize the cooling effects produced). Also, model problems involving different 
geometries (such as cylindrical or annulus) should be solved in a similar manner to the 
present parallel plate problem. A potential source of difficulty in solving problems using 
more complicated coordinate systems is the ability (or possible inability) of MAPLE to 
solve the final equation for the time-averaged temperature term. 
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